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Abstract 

This paper develops a two-dimensional structural framework for valuing credit default swaps 
and corporate bonds in the presence of default contagion. Modelling the values of related 
firms as correlated geometric Brownian motions with exponential default barriers, analytical 
formulae are obtained for both credit default swap spreads and corporate bond yields. The 
credit dependence structure is influenced by both a longer-term correlation structure as well 
as by the possibility of default contagion. In this way, the model is able to generate a 
diverse range of shapes for the term structure of credit spreads using realistic values for input 
parameters. 



1 Introduction 

Firms do not operate in isolation and company defaults are not independent. In reality a whole 
network of links exists between companies in related businesses, industries and markets and 
the impact of individual credit events can ripple through the market as a form of contagion. 
It is thus of fundamental importance when modelling credit, not only to understand the 
drivers of credit risk at an individual company, but also the dependence structure between 
related companies. Whether accounting for counterparty risk in the price of a single-name 
credit derivative, or considering credit risk in a portfolio context, an understanding of credit 
dependence is essential to accurate risk evaluation and pricing. 

Credit weakness, ratings downgrades and ultimately corporate default can occur in three main 
ways. Firstly, a company may be adversely affected for reasons specific to that company alone 
(e.g. poor financial management). Secondly, credit weakness may occur due to a factor or 
factors impacting multiple companies - whether in the form of a cyclical influence related 
to the economy, or a market- wide shock such as an earthquake or September 11 th . Finally, 
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companies are related through ties, some of which are real (e.g. a trade-creditor agreement), 
others of which are purely a matter of perception - for example the fear of accounting fraud. 
Credit dependence then occurs primarily through two mechanisms - either as a direct conse- 
quence of a common driving factor, or due to inter-company ties. The latter can be thought 
of as a form of contagion. 

The purpose of this paper is to examine the importance of considering a more complete de- 
pendence structure than is usually incorporated in credit models, one that better reflects 
reality, taking into account both a common driving influence and the possibility of idiosyn- 
cratic company links. To do so, whilst retaining some degree of analytical tractability, we 
consider a two-dimensional structural model with default as the first hitting time of an ex- 
ponential barrier. Firm values are modelled as correlated geometric Brownian motions and 
the default event is contagious. In this way, we are able to capture the two facets of the 
dependence structure. The correlation in firm values reflects a longer-term common driving 
influence on corporate strength whilst default contagion represents a direct link between the 
fortunes of both companies. As discussed in Section 3, this framework results in a model that 
is asymmetric with regard to default risk, a significant improvement on prior models. 

Structural models, whilst far from straightforward mathematically, are far more grounded 
in economic fundamentals than many other models and thus form a good starting point for 
a realistic description of credit dynamics. Giesecke (2004), Schonbucher (2003) and Lando 
(2004) provide a good introduction to structural models, their development since first intro- 
duced by Merton (1974), and their traditional place within credit modelling. Until recently, 
the vast majority of work on the structural model has focused on the case of a single firm. In- 
deed, two popular commercial packages, Moody's KMV and CreditGrades™, are motivated 
by the single-firm structural model 1 . Very little, however, has been published for multiple 
companies, with the market mainly focused on copulas or conditionally independent factor 
models 2 in the multivariate case. Two exceptions are the papers by Zhou (2001) and Hull and 
White (2001). Zhou (2001) calculates default correlations for two firms whose asset values are 
modelled as correlated Brownian motions. Hull and White (2001) extend this to framework 
to higher dimensions and proceed numerically in a discrete-time setting, proposing a method 
to calibrate piecewise constant default barriers to a term structure of default hazard rates. 
In contrast, we extend the framework used in Zhou (2001) to incorporate default contagion 
and derive analytical formulae for bond yields and CDS spreads. 

Copula methods, which allow the dependence structure of a portfolio to be considered inde- 
pendently from individual default times, are easy to implement and have rapidly become the 
market standard for modelling portfolios of credits. However, as basically static models able 
only to model expected defaults over a given time period, they fail to allow for suitable credit 
spread dynamics and tend to exhibit time instabilities. 3 Furthermore, the copula approach 
has no notion of default cause and effect as exists in a contagion mechanism. 

As problems have arisen with the widespread market use of copulas, multidimensional struc- 

1 Further details of these approaches can be found at www.moodyskmv.com and in Finger et al. (2002), 
respectively. 

2 For a good overview of the use of copulas in finance, see Cherubini et al. (2004); Schonbucher (2003) 
provides a useful summary and references for factor models. 

3 Mikosch (2006) provides a critical discussion of the widespread usage of copula methods; some countering 
arguments are given by Genest and Remillard (2006). 
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tural models are seeing renewed interest. Hull et al. (2005) price CDO tranches in a structural 
framework where assets are driven by a common factor. In this way, defaults are modelled in 
a dynamic setting, and firm value correlations can be time-dependent or stochastic, however 
the dependence structure stems purely from the correlated firm values and so is unable to 
account for any default causality or contagion. In another approach, Luciano and Schoutens 
(2005), Moosbrucker (2006) and Baxter (2006) assume that firm values are driven by Levy 
processes rather than geometric Brownian motions. By assuming that firm values are mod- 
elled as geometric Brownian motions time-changed by a common Gamma business time, firm 
values become Variance Gamma processes, allowing for a richer characterisation of spread dy- 
namics. Dependence is introduced through having a common stochastic time change, which 
may then be further broken down into systematic and idiosyncratic components, with a num- 
ber of different representations proposed by the various authors. The resultant dependence 
structure is more realistic and flexible but does not incorporate any form of default contagion. 

The paper is organised as follows. In Section 2 we outline the framework for the model, its 
assumptions and underlying results. We consider both the most general formulation and cases 
in which the formulae simplify. Formulae and results for corporate bond yields are provided 
in Section 3, whilst those for credit default swap spreads are in Section 4. We summarise and 
consider future extensions in Section 5. Mathematical details are in the appendix. 

2 The Model 

We consider two companies, firm values V{, i = 1,2. Each company issues equity and a 
single homogeneous class of debt, assumed to be a zero coupon bond, Ci(t,T), par value K%, 
maturity T. For simplicity we assume that both bonds have the same maturity date but the 
analysis is easily extendible to different maturity dates. 

For each company, firm value is assumed to follow a geometric Brownian motion, with default 
as the first time that the value of the firm hits a lower default barrier bi(t). As in Merton 
(1974) and Black and Cox (1976), we assume that a firm's value can be constructed from 
tradable securities and so in the risk-neutral pricing measure, for i = 1,2, 

dVi{t) = (r f - qi )Vidt + aiVidWiit) 

where the risk-free rate, rf, dividend yields, qi, and volatilities, Uj, are constants, Wi(t) are 
Brownian motions and cov(Vt / i(t), W2(t)) = pt for constant correlation p. 

We assume that each company has an exponential default barrier, reflecting the existence of 
debt covenants, and denote the default barrier for company i by 

Whilst similar to the barrier formulation in Black and Cox (1976) 4 , defining the barrier in 
this way provides us with additional flexibility, enabling us to change its slope. In particular, 
in the special case that the barrier growth rate is set equal to the drift in firm value, the 
model simplifies. As neither the barrier growth rate nor the drift in firm value is observable 

4 Black and Cox (1976) use a barrier of the form ui i K i e^ r f ( - T ^ t " > where < Wj < 1. The similarities and 
differences in this formulation compared to the one we use are discussed further in Section 3. 
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in practice, this is an attractive simplification. It is by no means necessary for our analysis, 
however, and we provide results for the most general case. 



Setting 



enables us to consider the simpler case of Brownian motion with drift and constant default 
barrier B { = In < 0. X;(0) = and 



where a* = ry - % - 7* - ±07. 
Defining the running minimum 



X i (t) = a i t + a i W i {t) 



X 4 (t)=mmX i (s), 

0<s<t 



and default time, Tj, as the first hitting time of the default barrier, 

Ti = mf{t: Xi(t) = Bi}, 

survival probability is then 

Pfa > s) = PQQ( S ) > Bi). 

The key result we use to value credit spreads is the joint survival probability density function 
and the resultant joint survival probability, P (t), 



P(t) = F(X 1 (t)>B 1 ,X 2 (t)>B 2 ) (1) 



where 



POO 

9n{0) = I re-^e^>I { ^(^)dr 

a±a 2 — pa 2 <Ji _ OL2d\ — pa\Oi 

ai ~ (1 - pi)o\o 2 ' ° 2 - (1 - p^ol 

b = — oi\ai — a2(i2 + 2 a \ a \ + P&i^a^ + -<r|a| 
fT~Z. — 2 

tan/? = - V - 9 , /3 G [0,vr] 

1 (B\ 2 P B X B 2 | j? 2 2 ^ 1/2 
V 1 - p 2 V °"i ^i^a ^2 

(T2-B1 - po\B2 
A{9) = a\(Ji sin(/3 — 6) + 0202 sin 6*, 



4 



and (^r 1 ) is a modified Bessel's function. P(t) represents the probability that neither 
company hits its default barrier by time t. Correlation between firm values is reflected by (3, 
with (3 = it corresponding to perfect correlation, (3 = tt/2 independence and (3 = perfect 
negative correlation. When the drifts in firm value and the default barrier are equal, a,i = 0, 
leading to ai = = b. The form of P(t) derives from the separation of variables in the 
solution of the Fokker-Planck equation governing the evolution of the survival probability 
density function. Full details of the calculation methodology can be found in the paper by 
Hua et al. (1998) in which double lookbacks are valued using the joint distributions for the 
maxima and minima of two correlated Brownian motions. Zhou (2001) uses the same result 
to calculate default correlation for two companies with correlated assets. In particular, Zhou 5 
primarily considers the special case that the barrier growth rate and drift in firm value are 
equal, 73 = rf — qi — , and shows that 



Pit) = ^e~^ V -sin(^) \n (m+1) (f) +I (if) 



„2\ 1 

• (2) 



This simplification, implying constant leverage, makes the joint survival probability consid- 
erably faster to evaluate. Zhou (2001) found that the assumption had minimal impact on 
default correlations for shorter maturities. Supporting this finding we see very little difference 
in survival probabilities as the barrier growth rate, ji, is changed for maturities up to five 
years. This is reflected in the fact that implied bond yields are not very sensitive to changes 
in the barrier growth rate, as illustrated in Figure 4. 6 

The joint survival probability calculation can also be simplified by selecting values of (3 for 
which the modified Bessel's function simplifies - for example (3 = | for integer k. Unfor- 
tunately, since p = — cos (3, the majority of these cases correspond to negative values of 
correlation and so are of only limited interest in our framework. 



3 Bond Yield Calculation 

The value of a corporate bond to a bondholder arises from two components - its value on 
maturity (should it mature) and its value in the event of default. We consider the t = value 
of a zero coupon bond issued by company one, par value Ki, maturing at time T, denoted 
Ci(0,T). The yield is then 

V*T) = -i.n(^). (3) 

We incorporate default contagion by assuming that company one defaults on its outstanding 
debt the first time that the value of either company reaches its default barrier. In this 

5 N.B. The definitions of ai and 0,2 in Zhou are the negative of the definitions in this paper. Our formula for 
g n (6) can be reconciled with Zhou's using double angle formulae and the fact that cos (3 = —p, sin (i = \J 1 — p 2 . 
(His angle a is our angle /3.) 

6 An extensive analysis of survival probability sensitivity to input parameters is contained in Haworth (2006) ; 
here we consider solely the impact of key parameter assumptions with regard to credit spreads, with results 
in Section 3. 



5 



way, there are two components to the default mechanism. If the value of firm one declines 
sufficiently, the company is forced into bankruptcy - this is due to the direct performance of 
the company itself and is exactly the framework used in normal first-passage structural models. 
The second way in which company one can default is if company two goes bankrupt, modelled 
as the time when the value of company two reaches its default barrier - default contagion. This 
would be the situation if company two was essential to the continuing operation of company 
one. For example if company two was the only purchaser of company one's products, as in 
the case of a small or regional auto-parts supplier to General Motors. 

Our specification of the contagion mechanism implicitly assumes that the full extent of its 
existence is not known to the market beforehand, 7 as is so often the case. Many links between 
companies are very opaque and certainly not fully disclosed. For example, the full nature of 
a bank's loan portfolio is rarely apparent, and there are many instances when banks have 
ended up over-exposed to one company or industry. The situations at Parmalat and LTCM 
in recent years are cases in point. Another example is the networks of business links that 
exist, for example in Italy, when extensive corporate cross-holdings are common. Few of these 
are publicly known. Too frequently, the existence or full extent of a company's relationships 
and the consequent investment risks only come to light when problems arise. The goal of our 
model is to take the first step in considering the spread implication of these types of links 
through the introduction of a contagion process that is simple enough to allow for analytical 
solutions. It is worth noting that company two need not default automatically if company 
one does. It can continue to operate regardless of the financial viability of company one with 
dependence solely through the asset correlation, p. As a result, the model is asymmetric with 
respect to default risk, in stark contrast with the majority of previous models incorporating 
a credit dependence structure. 8 

This framework is only really realistic for p > 0. It is highly unlikely that the bankruptcy of 
one company would lead to the immediate default of another, negatively correlated company. 
Whilst it is possible to contrive a theoretical example (e.g. a highly diversified company like 
General Electric might be key to one of its suppliers, but negatively correlated with it overall), 
economically it is rather improbable in practice and so we restrict ourselves to consideration 
of < p < 1. 

Similarly to Black and Cox (1976), we assume that 

Payment at maturity = min(a>iVi(T), K\) provided t\ > T, T2 > T 

where K\ is the par value of the bond, Tj denotes the default time of company i, and uj\ is a 
constant write down factor. This factor is the same as used later in the specification of the 
payment on default. It represents the fact that in the event of default or a restructuring, a 
portion of the defaulting company's value is lost to bondholders. This is commonly seen in 
practice, caused, for example, by restructuring costs or violation of the priority rule allocating 
claims in the event of default. 

7 We are grateful to one of our referees for highlighting this. 

8 The natural extension of this framework to larger portfolios of companies would be the type of situation 
considered in Jarrow and Yu (2001) in which 'primary' companies impact 'secondary' companies but not vice 
versa. A 'primary' company is likely to be larger with a greater market impact than a 'secondary' company. 
For example Microsoft or General Motors compared to a small, local IT or auto component manufacturer. 
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Black and Cox (1976) incorporate the write-down factor in the default barrier, but then the 
bondholders are always paid in full once the value of the firm at maturity is at least the face 
value of the bonds. This seems a little unrealistic since a company would be unable to pay 
out its entire value to bondholders. There would be liquidation costs and a company would 
not be able to raise its entire value in a refinancing. It therefore makes sense that the firm 
value at maturity must be equal to some Ki/r]i, where uj\ < rji < 1, for bondholders to be 
fully repaid. For simplicity, since fewer parameters are preferable, we assume that the firm 
repays bondholders in full for V\(T) > K\/u\. This would seem to be an improvement on 
the case where full repayment occurs for firm value at maturity of par or more. 

Changing to Xi(t) coordinates, the discounted maturity payment, DMP, can be written 



DMP 
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where 
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is the joint survival probability density function at maturity and 
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We assume that default occurs the first time that either company hits its default barrier 
and that in the event of default, the bondholder receives a percentage of discounted par 
value, uj 1 K l e~ r f ( - T - T \ This is the same payoff as used by Black and Cox (1976) and is highly 
attractive since the discounted default payment then becomes 

u 1 K 1 e- r f T {l-P(T)), (6) 

where the joint survival probability P{T) = P(2Ci(T) > Bi,X 2 (T) > B 2 ) is defined in (1). 

By construction, the default payment is worth less than discounted par, and so for consistency 
we just need to ensure that it is worth less than the value of the firm at default. Since company 
one must be worth at least as much as its default barrier, a sufficient condition is 

Ul < e (r/-m)(T-rO 

Adding together (5) and (6), and using (1) for P(T), yields can then be calculated using (3). 
In order to evaluate (5) and (6) we use the integral form of the modified Bessel's function, 

/rr \ 1 /twA 1 . (mr 2 \ f°° _im coshs _»™ 
I/H2L-) = — / e t LUS ^ cos — — \ dm sin — — / e * f> as. 

The infinite sums in (5) and (6) converge rapidly to zero, and so this substitution enables 
us to approximate solutions by a finite sum of three-dimensional integrals which we evaluate 
by numerical quadrature using a sparse grid (for further information regarding sparse grid 
methods, see Gerstner and Griebel (1998)). Figures 1-6 illustrate results and parameter 
sensitivities. As a measure of company strength we consider initial firm value, V$(0), divided 
by the initial level of the barrier, 6j(0). We denote this parameter by initial credit quality 
and spread sensitivity to it is as would be expected - as initial credit quality declines, spreads 
widen significantly. An indicative scaled distance to default, ^-log(-j^j), for our parameter 
values of (jj = 0.2 and initial credit quality of 2 is then 3.5. 

In all cases, we see that yields decline as correlation increases. Since default is less likely with 
increasing correlation, 9 the bond is less risky, bond-holders are not rewarded with such high 
returns, increasing the price and reducing the yield. 

Figures 1 and 2 illustrate the yield curve for two values of the write-down factor, u)±. Compar- 
ing them, we see that as u>\ decreases, the yield curve inverts. The extent to which payments 
are written down clearly has a large impact on yields, particularly at the short-end of the 
yield curve. The thick black line, labelled 'none' in the figures is the bond yield for a firm 
with the same parameter values, but operating in isolation. In other words, a single firm's 
yield when modelled in a first-passage framework, but with no exposure to another company 
through default contagion. As one would expect, yields are lower as such a firm's bonds are 
less risky. It is immediately apparent that allowing for default contagion has a significant 
impact on yields, especially for longer-dated bonds. 

Figure 3 shows the impact of varying u\ on a five-year bond. Of note, the case when uj\ = 1 
corresponds to no write-down on default and in effect the bond becomes risk-free, yielding, 

9 This is the case whether or not there is default contagion and can be seen from the fact that the probability 
of at least one of the companies defaulting in a given time period decreases as they become more correlated. 
If Ni denotes default by company i, for i = 1, 2, then P(JVi U N 2 ) = P(A?i) + P(iV 2 ) - P(iVi n N 2 ) and the final 
term increases with p. 
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Figure 1: Implied yield curve, uj\ = 0.7 



Figure 2: Implied yield curve, u\ = 0.5 





a x =a 2 = 0.2, K x = 100, r f = 0.05, 
9i = Q2 = 0, 71 = 72 = 0.03, initial credit quality = 2 



as we would expect, the risk-free rate of 5% regardless of correlation. uj\ acts in two ways - 
it lowers the payment in the event of default and it increases the value the company must be 
worth at maturity for bondholders to be repaid in full. 



Figure 3: Implied bond yield, varying oj\ 




0.4 0.6 
Correlation 



g { = 0.2, Ki = 100, r f = 0.05, qi = q 2 = 0, 
= 0.03, initial credit quality = 2, T=5 
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Figure 4: Implied bond yield, varying ji 




0.4 0.6 
Correlation 



a, 
1i 



0.2, K x = 100, r/ = 0.05, oji = 0.7, 
0, initial credit quality = 2, T=5 



Figure 4 considers the sensitivity of yields to the shape of the default barrier for a 5-year bond 
with initial credit quality of two and uj\ = 0.7. Changing the slope of the default barrier has 
minimal impact on yields - as the slope increases, default is less likely and yields decrease, but 
the impact is fairly small, particularly when considering the dependence on other parameters. 
7i has even less impact for shorter maturities, but becomes progressively more important as 
time to maturity increases. 

Finally, we consider the impact of varying the volatility of firm- value. Figure 5 shows how 
yields behave when the volatility of both firms is changed simultaneously. As expected, higher 
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Figure 5: Implied bond yield, varying <jj 



Figure 6: Implied bond yield, varying 02 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Correlation Correlation 

7i = 72 = 0.3, Ki = 100, t } = 0.05, ui = 0.7, 
qi = q2 = 0, initial credit quality = 2, T=5 years 



volatility leads to a higher likelihood of default and higher yields. In Figure 6 we assume that 
the volatility of firm one remains fixed at 0.2 yr -1 / 2 and we increase the volatility in firm 
two. Since we are considering the yield on firm one's bonds, the increasing riskiness of firm 
two impacts yields through the correlation between the two companies and the possibility of 
default contagion. As expected, the more volatile firm two is, the riskier firm one and the 
higher yielding its bonds. 

4 CDS Spread Calculations 

Using a similar approach to that in Section 3, we evaluate first and second-to-default credit 
default swap (CDS) spreads for a two-company basket. 

4.1 First-to-default CDS Basket 

We consider a basket of two related companies. The buyer of a first-to-default CDS on 
this underlying basket pays a premium, the CDS spread, for the life of the CDS - until 
maturity or the first default, whichever happens first. In the event of default by one of 
the underlying reference companies, the buyer receives a default payment and the contract 
terminates. Denoting the default time of company i by Tj, we write Tfi rs t for the time of the 
first default, 

r firs t = min{ri,r 2 } 
where, using the same notation as before, 

r i = mf{t:X i (t)=B i }. 

If bond recovery on default is R, and the protection buyer makes continuous spread payments, 
c, on a par value K, then the discounted spread payment (DSP) and discounted default 
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payment (DDP) on the first-to-default basket are 

DSP = cK [ e- r / s P(r first > a) ds 
Jo 

rT 



DDP = (1 - R)K [ e' r f s F(s < r first < a + ds) ds (7) 
Jo 

= (1 - R)K J -e- r f s —F( T{iTSt > s) ds 

= (1 - R)K^1 - e- r / T P(r first >T)-r f ^ e -^P(r first > s) da j. 

With P(a) defined as in equation (1), the market spread, Cfi rst , is therefore 

(1 - A) jl - e- r / T P(r first > T) - J T r /e ^ s P(r first > a) da j 

Cfirst = " " (8) 

/o e- r / s P(r first > a) da 

(1 - 12) jl - e~ r f T P(T) - r f e- r f s P(s)ds^ 
_____ 



since 



P(T firs t > s) = P(n > s,r 2 > s) = P(X x (a) > X 2 (a) > B 2 ) = P{s). 



4.2 Second-to-default CDS Basket 

A second-to-default CDS spread is evaluated in the same way. The purchaser of the swap 
receives a payment in the event that both companies default during the life of the swap, 
at which point the contract terminates. Denoting r seC ond as the time of the second default, 
exactly as for (8), the market spread, c secon( j is 



(1 - R)\l - e- r ' T P(r sccond > T) - / T r /e -^P(r sccond > a) da j 

Csecond = ~7r ', [ ' (^) 

/o e- r f s ¥(r second > a) da 



where 

P(r SC cond > s) = P(n > a) + P(r 2 > a) - P(n > a, r 2 > a). (10) 
4.3 CDS Basket Results 

In Figures 7 - 12, we consider the impact of correlation on first and second-to-default CDS 
spreads for different parameter values. Numerical evaluation is done by numerical quadrature 
on a sparse grid as before. 
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Figure 7: First-to-default CDS, varying T 



Figure 8: Second-to-default CDS, varying T 




0-1=0-2 = 0.2, Ki = 100, rj = 0.05, q x = q 2 = 0, 
7i = 72 = 0.03, initial credit quality = 2, R = 0.5 



In all cases, with increasing correlation between the two reference entities, first-to-default CDS 
spreads decrease, whilst second-to-default CDS spreads increase. This is because the proba- 
bility of at least one company defaulting in a given period is higher for negative correlations, 
whilst the probability of both defaulting is greater for positive correlations. 

Figures 7 and 8 show spreads for first and second-to-default CDS baskets with maturities 
of up to 5 years. Initial credit quality is 2 (i.e, as before, firm value is initially twice the 
level of the barrier). Spreads are greater for longer-maturity swaps and, as we would expect, 
first-to-default spreads are everywhere greater than second-to-default spreads. 

Figure 9: First-to-default CDS, varying R Figure 10: Second-to-default CDS, varying R 




Correlation Correlation 



0-1=0-2 = 0.2, K x = 100, Tf = 0.05, qi=q 2 = 0, 
7i = 72 = 0.03, initial credit quality = 2, T = 5 

Figures 9 and 10 illustrate the extent to which CDS spreads depend on our recovery rate as- 
sumption. As would be expected, moving from a 30% recovery rate to a 70% recovery rate has 
a large impact. However, the overall form of spreads and their variation with changing cor- 
relation is the same. In general, taking R=50% is representative of the levels seen in practice 
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(see, for example, Bakshi et al. (2006)) and is in line with that used in the CreditGrades 
approach to modelling credit as described in Finger et al. (2002). An easy extension would 
be to set the recovery rate of the reference entities equal to discounted par value. 




Figures 11 and 12 show that firm volatility has a considerable impact on spreads. As the 
reference entities become more volatile, both credit default swaps become much more risky 
and spreads increase significantly. Of all the parameters, firm volatility has by far the largest 
impact on CDS spreads. 

Through implementation of the analytical formula for the joint survival probability function, 
(1), we are therefore able to illustrate the sensitivity of first and second-to-default CDS 
spreads to input parameter assumptions straightforwardly. The importance of the degree of 
correlation between reference entities is clearly evident, with spreads significantly different 
when the basket is highly correlated than when it is well diversified. 

Since by definition default contagion comes into play following the first default, its existence 
has no impact on first-to-default spreads. It is, however, important in the consideration of 
second-to-default spreads, and introducing the default contagion mechanism considered in 
Section 3 clearly has a dramatic impact on second-to-default spreads, in this case making 
them equal to first-to-default spreads. This situation is clearly extreme, as a CDS is highly 
unlikely to be constructed on such an undiversified basket (unless of course the existence of 
the link between firms was completely unknown prior to the first default), however it serves 
to highlight the importance of considering the dependence relationship between firms. 



4.4 CDS with Counterparty Risk 

Consider now a single-name CDS, face- value K, maturity T, on reference company one bought 
from a counterparty company two. The purchaser of the CDS makes spread payments for 
the life of the CDS - until either the reference company or the counterparty defaults. If the 
reference entity defaults during the life of the CDS and before the counterparty, the purchaser 
receives a default payment. If, however, the counterparty defaults first, they receive nothing, 
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irrespective of whether or not the reference company later defaults. Denoting the default 
time of company i by Tj, if bond recovery on default is R and the purchaser of protection 
makes continuous spread payments, c, for the life of the CDS, then using the same notation 
as before, the protection buyer 

• makes spread payments for t < minjri, t 2 , T}, 

• receives a default payment if t\ < min{r2,T}, 

• receives nothing if t 2 < min{ri,T}. 

In other words, discounted spread and default payments are 



Spread = ck[ e" r ^ s P(ri > s, r 2 > s) ds 
Jo 

= cK [ e~ r f s ¥(X 1 (s)>B 1 ,X 2 (s)>B 2 )ds 
Jo 

Default = (1-R)K^J e^' s P(s < n < s + ds,T 2 > s) ds 
Considering the default payment for a first-to-default CDS, equation (7), 

DDP = (1 - R)K f e~ r f s F{s < r first < s + ds) ds (11) 
Jo 

= (1-R)kJ^ e~ r -f s |p(s < n < s + ds,r 2 > s) +P(n > s,s < r 2 < s + ds)| ds, 

we see that the default payment for a CDS with counterparty, added to its image when the 
identity of the reference entity and the counterparty are swapped, gives the value of the first- 
to-default swap payment. A similar identity holds for the second-to-default swap payment. 

In the case of a homogeneous portfolio, when both reference entities have the same parameters, 
(11) can be used to calculate the value of a CDS spread with counterparty risk for all values 
of correlation, p. More generally, in the asymmetric case a more complicated evaluation must 
be done, details of which will be considered elsewhere. 

5 Conclusion 

Structural models are increasingly the focus of industry attention in the multi-firm setting, 
but there has been limited work on the general first passage framework, and little academic 
coverage. This paper makes a first and novel, albeit by necessity simplistic, contribution 
to this field through the introduction of a contagion mechanism in a two-dimensional first 
passage model. Working with a Black and Cox (1976) type structural framework, we have 
built on the work by Zhou (2001) to derive analytical formulae for both bond yields and CDS 
spreads. We have modified the default barrier to better reflect reality and have incorporated 



14 



default contagion within the structural framework for the first time. The result is a credit 
model that is asymmetric with respect to default risk and which has a dependence structure 
based on both long-term asset correlation and default contagion. 

Results illustrate that the sensitivity of yields to input parameters is as expected, and clearly 
demonstrate the importance of credit correlation. Our model enables us to generate corporate 
bond yields and CDS spreads across the full range of parameter values in a two-dimensional 
first passage framework in full generality. This has not been done before. Previous results 
using related analysis in Zhou (2001) (default correlations) and Hua et al. (1998) (double 
lookbacks) have concentrated on cases in which the framework simplifies and have been lim- 
ited to a few parameter values in these cases. For the first time, we are therefore able to 
fully consider spread sensitivity to model and parameter assumptions in the two-dimensional 
structural setting. 

Our specification of default contagion is clearly not very realistic - default by one company 
very rarely leads to direct default by another, although it is possible. More likely, the im- 
pact of a corporate bankruptcy causes a ripple of credit weakness through the market as 
related companies are impacted. 10 Nonetheless, the importance of taking into account credit 
interactions is, once again, clearly highlighted. 

Dependence modelling is most critical in the analysis and pricing of large basket credit deriva- 
tives, such as k ift -to-default credit default swap baskets and CDO tranches. These require 
the framework to be extended to considerably more than two dimensions. This is an area of 
current research interest and is not straightforward since analytical formulae are no longer 
possible and numerical solutions become highly problematic with increasing dimension. As 
intimated in Section 3, it would be attractive to incorporate a network of asymmetric depen- 
dences within a portfolio, enabling the impact of a credit event at one company to cause a 
ripple of contagion through other, related, parties. 

A Derivation of Maturity Payment 

For ease of notation, we denote the joint survival probability transition density 

d 2 

p{x 1 ,x 2 ,t) = dxidx2 F ( X i(t) < x u X 2 (t) < x 2 ,X 1 (t) > B 1 ,X 2 (t) > B 2 ) 
Proposition A.l For general A and B, 

roo rB 

/ / e eXl p(xi,x 2 ,t)dx 1 dx 2 

= 1. e (a 1+ e)B 1+ a.B 2+ M £ e -rg/2t gin f f ^ ( ™#\ ^ ^ 

pt n=1 V p J Jo V p ) 

1() We address this numerically in a subsequent paper, Haworth and Reisinger (2006), however analytical 
solutions are no longer possible. 
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where 



9n(0) 

tan 6*o 

d A 



d B 



dA 



re -r 2 /2t e [A(tf)+ e <n B in03-(»)]r J( ^ ) dr 

1 (B\ 2pB 1 B 2 Bl\ 1/2 



tJI - p 2 V °\ a i a 2 
a 1 B 2 y/l - p 2 



+ 



4 



a 2 Bi - pcr 1 B 2 ' 
A — Bi 



O € [0,/3] 



0"! 



\J\ — p 2 cos 6* + p sin 
B - B x 



d B 

A(9) = aia\ sin(f3 — 9) + a 2 a 2 sin 9 , 



\J\ — p 2 cos 9 + p sin 



Proof 

Using notation from Section 2, by Hua et al. (1998) 



p(xi, X 2 , t) = 7 == 



_C r 2 +r 2y 2t . I nir9o\ . ( nirv\ T /rr \ 



Changing variables, 



xi = Si + -y/(l — p 2 )o~\r cos 6* + po^r sin I 
x 2 = B 2 + a 2 rsin9, 



the Jacobian for the transformation is y/ (1 — p 2 )ra\a 2 , and 



roo i-B 

I I e tXl p{xi,x 2 ,t)dxi&x 2 

JB 2 J A 



^Ja l +e)B 1 +a 2 B 2 +USr p -rl/2t • f n7r6l 



sin | ^ ) g n {0)d9 



where 



since 



n (0) = ^ re -r 2 /2 ie [A W+ e CT1 si„( / 3- e) ]r /( _ ) ^ 



cos (3 = —p & sin /3 = y 1 — p 2 . 
Writing X = and y = then from (12), 



X = r 



a/(1 — p 2 ) cos + p sin 9 



Y = r sin ( 
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Since 



A - B l B-B 1 

0< < X< 

CJl (7l 

Y < oo 



it follows that 



and letting 



tan 



X/Y-p> 



0e[0,(3] 



d A 



A-B x 



\J\ — p 2 cos 9 + p sin 9 



B-Ba 



\/l — p 2 cos 9 + p sin 9 



we have d A < r < ds from (13). The result follows. 



Derivation of Maturity Payment 

From Section 3, the discounted maturity payment, DMP, is: 

/'OO rco 

DMP = e _r ' T I I K 1 p(x 1 ,X2,T)dx 1 dx 2 

JB 2 Jd 

poo rd 

+ e~ r i T / u 1 V 1 (0)e x ^ T p(x 1 ,x 2 ,T)dx 1 dx 2 



□ 



(14) 
(15) 



Using Proposition (A.l) with e = 0, A = d and B = oo for line (14) and e = 1, A = B\ and 
B = d for line (15), the payment on maturity becomes: 



DMP = H l( T)± ^(^^sin^rfWdO 

+ H 2 (T)± sin (=5*) jf rtWd* 



17 



where, 



Hi(T) = 




e a 1 B 1 +a 2 B 2 +bT g-rg/: 



2u 1 V 1 (0)e^- r ri T 
/3T 



{a 1 +l)B 1 +a 2 B 2 +bT e -rl/2T 



r 

= / 

Jd 

am = I 

Jo 



oc 



/2TA(8)rr frr \ 



re 



dr 



d*(d) 
,d*(0) 



-r 2 /2T e [A(0)+«ri sin(/3-e)]r r 



/rr \ 

(^) 



re 



dr 



d- B 



i 



<7i a/1 — p 2 cos + p sin 
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